New Universality Classes for 
Two-Dimensional cr-Models 



Sergio Caracciolo 
Scuola Normale Superiore and INFN - Sezione di Pisa 
Piazza dei Cavalieri 
Pisa 56100, ITALIA 
Internet: CARACCIOOUXISIS . SIS . IT 
Bitnet: CARACCIOOIPISISVA.BITIET 
Hepnet/Decnet: 39198 :: CARACCIOLO 

Robert G. Edwards 
Supercomputer Computations Research Institute 
Florida State University 
Tallahassee, PL 32306 USA 
Internet : EDWARDSOMAILER . SCRI . FSU . EDU 

Andrea Pelissetto* 
Alan D. Sokal 
Department of Physics 
New York University 
4 Washington Place 
New York, NY 10003 USA 
Internet: PELISSET@MAFALDA.PHYSICS.IYU.EDU, S0KAL@ACF4.IYU.EDU 



July 31, 1993 



PACS number(s): ll.lO.-z, 05.50.+q, 11. 15. Ha, 64.60. Fr 



*Address until August 31, 1994. Permanent address: Dipartimento di Fisica and 
INFN - Sezione di Pisa, Universita degli Studi di Pisa, Pisa 56100, ITALIA. Inter- 
net: PELISSET@SOTTHPI1.DIFI.OTIPI.IT; Bitnet: PELISSET@IPISISVA.BITIET; Hepnet/Decnet: 
39198: : PELISSETTO. 



Abstract 



We argue that the two-dimensional 0(iV)-invariant lattice cr-model with 
mixed isovector/isotensor action has a one-parameter family of nontrivial con- 
tinuum limits, only one of which is the continuum ci-model constructed by 
conventional perturbation theory. We test the proposed scenario with a high- 
precision Monte Carlo simulation for N = 3,4 on lattices up to 512 X 512, 
using a Wolff-type embedding algorithm. The finite-size-scaling data confirm 
the existence of the predicted new family of continuum limits. In particular, 
the RP^~^ and iV-vector models do not lie in the same universality class. 



Two-dimensional nonlinear cr-models are of direct interest in condensed-matter 
physics, and they are important "toy models" in elementary-particle physics by virtue 
of their close similarity to four-dimensional gauge theories. In this Letter we summa- 
rize analytical and Monte Carlo evidence [1] which calls into question some aspects of 
the conventional wisdom regarding two-dimensional cr-models and their universality 
classes. 

The conventional wisdom [2] holds that two-dimensional lattice cr-models with a 
non-Abelian symmetry group are asymptotically free — i.e. the theory has a mass gap 
for all /3 < oo — and that the asymptotic behavior as /3 — )• oo can be deduced from the 
lattice Hamiltonian simply by noting its symmetries and computing its perturbative 
expansion around an ordered state. The first part of this conventional wisdom has 
been questioned by Patrascioiu and Seller [3]. Here we do not question the first part, 
but rather the second part: we argue that there are new universality classes which 
cannot be seen in conventional perturbation theory. 

Consider the lattice cr-model taking values in the unit sphere S^~'^ C IR^, with 
Hamiltonian (= Euclidean action) 



(the so-called mixed isovector/isotensor model). According to the conventional wis- 
dom, this is simply a "variant action" for the usual (9(A^)-invariant S^~'^ cr-model, 
and all paths in the /3x)-plane satisfying /3 = /3y + /3x — ^ +oo are claimed to 
give rise to the same continuum-limit theory. Here we shall argue that this is not so: 
we claim that in fact there is a one-parameter family of continuum-limit theories, la- 
belled by a parameter \rem only one of which {\ren = +oo) is the continuum cr-model 
constructed by conventional perturbation theory. 

To see this, recall first that ai fij = -\-oo the spins are forced to point along a 
single axis; thus, the isotensor correlation functions are identically 1, and the isovector 
correlation functions are those of the Ising model at inverse temperature fiy- In 
particular, there is a critical point at j3j = +oo, /3y = /Seising = |log(l + a/2) = 





0.440686 . . . . 
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Figure 1: Proposed qualitative form of the "curves of constant physics" in the 
/3T)-plane. 



Assuming that there is no critical point in the finite /3x)-plane, let us inquire 
about the shape of the contours of fixed isotensor/isovector mass ratio niT jniv- More 
generally, for any dimensionless combination O of long-distance observables, let us 
inquire about the shape of the contours of fixed 0; asymptotically in the continuum 
limit these contours are expected to be independent of the observable 0, and we call 
them "curves of constant physics". 

We shall now argue that the curves of constant physics have the qualitative form 
shown in Figure 1. For concreteness, let O be (for now) the mass ratio niT/mv- 
For /3y = (the RP^~^ model), we have my = 00 by Z2 gauge invariance, hence 
niT/mv = 0. More generally, along any path satisfying fij ^ 00 and fiy < Pcjsing — 
we have mx — ^ while 1/my remains bounded above by its Ising value (by a rigorous 
correlation inequality [1]), so that niT jniv — ^ 0. On the other hand, it is heuristically 
reasonable that along any path satisfying /3 = /3y + /3x — ^ 00 and fiy > Pcjsing + 
the continuum limit should be exactly that of the A^- vector model, and in particular 
should have niT jniv — ^ 2. This is suggested qualitatively by the fact that at fij = +00 
and fiv > Pcjsing^ the isovector sector is spontaneously magnetized and has a finite 
correlation length; therefore, we expect that near such a point (at fij finite but large) 
the only possible critical fiuctuations are spin waves around this ordered state, which 
are identical to those in the pure A^-vector model. A more quantitative support for 
this view is provided by the expansion in powers of 1//3t around the Ising model 
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[4], which for /3y fixed > ficjsing is essentially identical to that of the pure A^- vector 
model, modulo a finite renormalization of fi. Putting together the two halves of 
the argument, we conclude that all the curves of constant physics must approach the 
fij = +00 axis precisely at fiy = Pcjsing- 

It follows from this picture that there exists (at least) a one-parameter family of 
continuum-limit theories, which can be parametrized by any suitable observable O. 
At one extreme is the RP^~^ continuum theory, which is obtained along any path 
satisfying fij ^ 00 with fiy < Pcjsing — £• At the other extreme is the A^-vector 
continuum theory, which is obtained along any path satisfying /3 = /3y + /3x — ^ 00 
with fiv > Pcjsing + £• Finally, there is a one-parameter family of intermediate 
theories, which are obtained along very special paths fij — )• 00, fiy Pcjsing that are 
asymptotically tangent (to sufficiently high order) to one of the curves of constant O. 

A renormalization-group argument suggests that the curves of constant O have 
asymptotically the form 

fiv = fifh) + Kengifh) , (2) 

where f{fiT) is a nonuniversal function arising from analytic corrections to scaling 
("nonlinear scaling fields"), g{fiT) is a simple universal function determined by the 
ratio of RG eigenvalues for the two relevant couplings at the Ising fixed point (these 
couplings are essentially fiy — Pcjsing and I//3x), and Xj-en is a parameter labelling the 
various curves. Then any observable is a universal function of Xj-en (modulo one 
or two nonuniversal scale factors), and we can label the possible continuum limits by 
their values of Xj-en'- 

(a) Xren = — oo corrcspouds to the RP^~^ continuum theory: it is controlled by 
the ^^RP^~^ fixed point", which is Z2-gauge-invariant and has one marginally 
unstable direction (basically Ij fij). 

(b) Xren = +00 corrcspouds to the conventional continuum cr-model: it is controlled 
by the 'W- vector fixed point", which has one marginally unstable direction 
(basically 1//3). 

(c) —00 < Xren < +00 corrcspouds to the one-parameter family of intermediate 
theories: they are controlled by the "critical Ising fixed point", which has one 
unstable direction with exponent v = 1 (basically fiv — Pcjsing) and one direction 
which is either weakly unstable or marginally unstable (basically Ij fij). 

This picture is confirmed in the Migdal-Kadanoff approximation, which predicts 
g{I^T) = e~™"'^*(^)^'^^ [I]. It is also confirmed by the exact solution in the limit 
A^ — )• 00 with fiv and fij = fij I N fixed [1, 4]. A more exact analysis can presumably 
be obtained by studying the renormalization properties of the expansion in powers of 
I//3t [4] around the Ising critical point fiy = Pcjsing] we have not yet succeeded in 
doing this, but we do note that for fiy = Pcjsing there are new infrared divergences 
beyond those in the A^- vector model, confirming our belief that new continuum the- 
ories — different from those of the RP^~^ and A^- vector models — can be obtained 
by suitable limits /3x — ^ 00, /3y — ^ Pcjsing- 
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We remark that the curve niT/mv = 2, which is the threshold for the generation 
of isotensor bound states, is one of the "curves of constant physics", corresponding 
to some particular finite value X^en = For X^en > the mass ratio niT/mv 
is identically 2, but the scattering matrix of the isovector particles is presumably a 
nontrivial function of X^en- 

Our family of continuum theories has a very peculiar feature: at the two extreme 
ends — the RP^~^ and A^-vector models — the isotensor correlation functions are 
identical, after a suitable change of length scale and redefinition of the coupling, at 
all orders in perturbation theory in Nevertheless, at the nonperturbative level 

it seems clear that the isotensor correlation functions cannot be identical in the two 
theories: the isotensor two-point function of the A^-vector theory presumably has a 
cut beginning at mass mx = 2my (i.e. Euclidean momentum p = zmx), corresponding 
to the scattering states of two isovector particles; while that of the RP^~^ theory 
presumably has a pole at mx, corresponding to a true isotensor particle ("bound state 
of two confined isovector constituents"). Apparently, nonperturbative corrections of 
order ^-'^'^f^/iJ^-'^) — which in the continuum theory correspond to corrections of order 
{p/^Ms)~'^ X logs, and are thus small at large momenta — contribute so strongly at 
low momenta as to change completely the particle structure of the theory. If this is 
in fact the case, then the perturbative A-parameter formula [5, 6] gives the change of 
length scale needed to match the correlation functions of the RP^~^ and A^-vector 
models at high momenta (i.e. much less than the ultraviolet cutoff but much greater 
than mx), but it says nothing about the behavior at small momenta and in particular 
about the particle masses. In other words, the conventional prediction [7] for the 
nonperturbative constant rnj /^ms iii the RP^~^ theory — namely that it takes the 
same value as in the A^- vector theory, viz. twice my/Aj^^ (which is now known exactly 
[8]) — is apparently invalid. 

The same argument shows that the intermediate theories have a nontrivial de- 
pendence on Xrem even in the isotensor sector, at least for X^en ^ Indeed, in 
any of the theories with an isotensor bound state (0 < niT/mv < 2), the isotensor 
correlation function has presumably a pole at my and a cut at 2my; so two theories 
with different mass ratios niT/mv cannot possibly be equivalent modulo a length 
rescaling. Almost certainly this nontrivial dependence on X^en continues even in the 
regime X^en ^ 

We have carried out extensive Monte Carlo runs on the mixed isovector /isotensor 
model (1) for A^ = 3 and N = on L X L lattices up to Z = 512 (resp. L = 128), 
using a Wolff-type embedding algorithm [9] . For A^ = 3 we have traced out the entire 
band-shaped region 5 < Cmax = ^^^{CvtCt) ^ 100-200: this involved runs at 558 
different triplets {fiy^fiT^L), using a total CPU time of about 7 years on high-speed 
RISC workstations. We measured the susceptibilities xy, xt and the second-moment 
correlation lengths iv\ij\ among other observables. Our lattice sizes typically 
satisfy 3 ^ L/^max ^ 10, and we have carried out an extensive finite-size-scaling 
analysis. Details of the static and dynamic behavior can be found in [1] and [10], 
respectively; here we summarize the highlights of the former. 

We constructed finite-size-scaling (FSS) curves for the observables O = , , 
Xy, Xt by considering all data points {fiy^fiT^L) for which we have data also at 
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/3t, 2iv), and plotting 0{2L) /0{L) versus ^max{L) / L. All points corresponding to 
the same continuum-limit theory should then lie on a single curve (modulo corrections 
to scaling); if our claims about the universality classes are correct, these universality 

(2) (2) 

classes can be parametrized by (for example) the ratio Cv /Ct ■ We thus predict 



0{2L) 



Fo{UAL)/L; ei'H^)/#(X)) + 0{L--) . (3) 



(2) (2) 

For the isovector (resp. isotensor) observables, we have divided {Cv '/Ct ') 

-space into 12 

(resp. 9) "slices" in which the scaling function Fq is found empirically to be constant 
within error bars, and in each slice we have chosen a smooth fitting function. The 
full set of curves is shown in [1]. The quality of the fits is amazingly good, and the 
curves vary smoothly from one slice to the next. More strikingly, the shape and even 
the sign of the finite-size corrections varies radically as a function of i^y VCt^- 

These results provide strong evidence for the scenario proposed above, in two 
ways: First, the fact that all the data points with the same value of i^y VCt ^ f^-U quite 
accurately onto a single curve — and this for all four observables — is strong evidence 
that (a) we are in the scaling region, and (b) all the theories with the same value of 

(2) (2) 

Cv /Ct belong to the same universality class. Conversely, the fact that the curves 

(2) (2) 

for different values of Cv /Ct ''^^^ same is strong evidence that theories with 
different values of i^y VCt^ belong to different universality classes. 

For the universality class of the A^- vector theory, our FSS curves for and xv 
are in close agreement with those predicted by Flyvbjerg and Larsen [11] from the 
1/A^ expansion through order 1/N'^. 

The contours of constant i^y VCt^ shown in Figure 2. These curves confirm 
beautifully the scenario of Figure 1: it is perfectly believable that all these curves 



(2) (2) 

are heading for the Ising critical point /3x = +oo, /3y 0.44. The ratio Cv /Ct 
essentially constant at its asymptotic value 3.44 in the entire region to the right of 
the last curve shown; in this area the theory is fully 'W- vector-like" . 

For the pure RP^~^ models with N = 3,4,6,8, our simulations up to correlation 
length ~ 225, 25, 10, 10 do not reveal asymptotic scaling. Therefore, we are unable 
to make any reliable statement about the nonperturbative constant mx/Aj^^. 

Our data for the RP^ (resp. RP^) model are also consistent with a different 
scenario [12], in which there is a critical point at /3x ~ 5.73 (resp. 6.96), only 3% 
(resp. 5%) beyond our last run. But for theoretical reasons we believe that such 
a critical point is unlikely: if it existed, then highly plausible (though unproven) 
correlation inequalities [1] would imply the existence of a critical curve continuing to 
the A^- vector axis (/3x = 0), in contradiction with the believed (though unproven) 
asymptotic freedom of these models [13]. Moreover, the critical-point scenario works 
less well when we try to trace the presumed critical curve into the (/3y, /3x)-plane 
and to estimate the critical exponents. We think it more plausible that the apparent 
critical point arises from a rapid rise of in a crossover region. 

The theoretical scenario proposed in this Letter applies also to other cr-models, for 
example the SU{N) chiral model with mixed fundamental/ adjoint action [14]. For 
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Contours of Constant ^v/^t 



5 ■ 
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Figure 2: Contour lines of constant ratio Cv j Ct = 0-3, 0.5, 0.75, 1.0, 1.25, 1.5, 2.0, 
2.5, 3.0, 3.2 in the /3x)-plane: expanded view of the region 0.4 < /3y < 1.0. Dots 
indicate our runs. 
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< 4 we predict the same scenario as in Figure 1, with new universahty classes con- 
trolled by the critical Zjv fixed point. For > 4 the situation gets more complicated, 
as the Zjv models are then believed to have an intermediate massless (spin-wave) 
phase [15] (this is proven for large N [16]). We conjecture that in this case the curves 
of constant physics end up on the /3x = +oo axis precisely in the intermediate-phase 
interval, with m^/mj? = running to the first (Kosterlitz-Thouless) transition point 
and niAlrnp = 2 running to the second (freezing transition). Since the difference be- 
tween SU{N) and SU{N) / is invisible in conventional perturbation theory, these 
scenarios are in confiict with the conventional wisdom. 

We predict a similar scenario for the complex mixed isovector/isotensor model, 
which interpolates between the 2A^-vector and CP ~^ models. In this case we con- 
jecture that the curves of constant physics end up on the /3x = +oo axis between 
/3y = /3c,XY and /3y = -|-oo; this picture is consistent with the large-A^ calcu- 
lations of Campostrini and Rossi [17]. Here, however, the manifolds S'^^~^ and 
CP ~^ ~ S'^^~^ /U{1) are different at the perturbative level; so our prediction, 
though apparently new, does not contradict the conventional wisdom. 

Similar scenarios can be expected for lattice gauge theories with mixed fundamen- 
tal/adjoint action [18] whenever the gauge theory based on the center of the original 
gauge group has a second-order (critical) transition. 

The results of this Letter show that the complete information about a continuum 
quantum field theory cannot be deduced solely from its formal continuum action. 
Rather, it is necessary to impose an ultraviolet cutoff (e.g. choose a lattice action), 
search the phase diagram of the cutoff theory for critical points, and study the possible 
continuum limits. We suspect that there are many hitherto-unexpected continuum 
quantum field theories waiting to be discovered. 
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